A generalized Caroli formula for transmission coefficient with lead-lead coupling 
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We present a generalized transmission coefficient formula for the lead-junction-lead system, in 
which interaction between the leads has been taken into account. Based on it the Caroli formula 
could be easily recovered and a transmission coefficient formula for interface problem in the ballistic 
system can be obtained. The condition of validity for the formula is carefully explored. We mainly 
focus on heat transport. However, the corresponding electrical transport could be similarly dealt 
with. Also, an illustrative example is given to clarify the precise meaning of the quantities used 
in the formula, such as the concept of the reduced interacting matrix in different situations. In 
addition, an explicit transmission coefficient formula for a general one-dimensional interface setup 
is obtained based on the derived interface formula. 

PACS numbers: 05.70.Ln, 44.10.+i, 63.22.-m 



I. INTRODUCTION 

In recent years there has been a huge increase in re- 
search and development of nanoscale science and tech- 
nology, with the study of energy and electron transport 
playing important role. Focusing on thermal transport, 
Landauer-like results for steady-state heat flow have been 
proposed earlier [l|, Subsequently, based on quan- 
tum Langevin equation approach, many authors success- 
fully obtained a Landauer-type expression . Alterna- 
tively, noncquilibrium Green's function (NEGF) method 
has been introduced to investigate mesoscopic thermal 
transport, which is particularly suited for use with bal- 
listic thermal transport and readily allows the incorpora- 
tion of nonlinear interactions |6|-|8| . Generally speaking, 
in the lead-junction-lead system, steady-state heat cur- 
rent of ballistic thermal transport flowing from left lead 
to right lead has been described by the Landauer-like 
formula, which was derived first for electrical current, as 



2^ 



(1) 



where f^L.B.} ~ {cxp {fi^/kBT^L,R}) ~ 1} is the Bose- 
Einstein distribution for phonons, and T [lo] is known as 
the transmission coefficient. Based on noncquilibrium 
Green's function method, T [uj] can be calculated through 
the Caroli formula in terms of the Green's functions of 
the junction and the self-energies of the leads, 

T[oj]=Ti-{G''TrG^Tl), (2) 
where G""'" is the Green's function of the junction, and 
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where the self-energy terms S^'^ are due to the semi- 
infinite leads on the left, L, and on the right, R, respec- 
tively. The superscript r and a denote the retarded and 
advanced, respectively, both for the self-energies as well 
as for the Green's functions in the formula. The specific 
form © was given from NEGF formalism by Meir and 
Wingreen for electronic case and later by Yamamoto 
and Watanabe for phonon transport [l^ , while Garoli et 
al. first obtained a formula for the electronic transport 
in a slightly more restricted case Also, Mingo et al. 
have derived a similar expression for transmission coeffi- 
cient using "atomistic Green's function" method [lH, [l^ . 
Very recently. Das and Dhar 14| derive the Landauer- 
like expression from plane wave picture using Lippmann- 
Schwinger scattering approach. 

The Landauer-like formula describes the situation in 
which the junction is small enough compared to the co- 
herent length of the waves so that it could be treated as 
elastic scattering where the energy is conserved. Further- 
more, it has been assumed that the two leads are decou- 
pled which physically means there is no direct tunneling 
between the two leads. Through modern nanoscale tech- 
nology, small junction is easily realized such as in certain 
nanoscale systems, for instance, a single molecule or, in 
general, a small cluster of atoms between two bulk elec- 
trodes. In that case, the electrode surfaces of the bulk 
conductors may be separated by just a few angstroms 
so that some finite electronic coupling between the two 
surfaces is inevitable taking into account the long-range 
interaction. In order to solve this problem, Di Ventra 
suggested that [l^ we can choose our "sample" region 
(junction) to extend several atomic layers inside the bulk 
electrodes where screening is essentially complete so that 
the above coupling could be negligible. It turns out to be 
correct using this trick to avoid the interaction between 
the two leads, which will be verified in a simple example 
at the end of the paper, even though we, to some limited 
extent, modify the initial condition necessary to derive 
Landauer-like formula in NEGF formalism and reparti- 
tion the total Hamiltonian. However, this procedure or 
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trick could not be always done due to some topologi- 
cal reason such as studying heat current in Rubin model 
[l6l | in which the other end of the two scmi-infinitc leads 
is connected (a ring problem). Actually this somewhat 
trivial example is not so artificial since it is equivalent to 
using periodic boundary condition in Rubin model. Fur- 
thermore, the modification of the initial product state 
will certainly affect the behavior of the transient heat 
current. If we want to study the transient and steady 
heat current |17| in a unified way, the repartitioning pro- 
cedure which changes the model is not acceptable. So in 
this work we will try to derive a compact formula applica- 
ble to this general model including lead-lead interaction 
for steady-state heat current according to NEGF formal- 
ism, and correspondingly obtained a Caroli-like formula 
for transmission coefficient. Furthermore, an interface 
transmission coefficient formula in the NEGF formalism 
will be given as a special case of the general Caroli-like 
formula. Also, the standard Caroli formula follows as a 
one-line proof. 

The paper is organized into two main sections. In Sec- 
tion ini we develop our formalism to derive our general- 
ized expression for the steady current directly taking cou- 
pling between leads into account. Based on this general 
formula, we will recover the Caroli formula and derive a 
computationally efficient interface formula in III Dl Then 
we apply this formalism to an illustrative model system 
in Section IIIII and show the results of numerical calcula- 
tions. Also, we will apply the interface formula obtained 
in Section In] to derive an explicit expression for transmis- 
sion coefficient in Section ITVl Finally we conclude with a 
short discussion in Section V. 



II. FORMALISM 

A. Model system 

As was mentioned previously, we will consider the lead- 
junction-lead model initially prepared in product state 

P(.to) = T.(e-^.«.) ® Tr(e-^C^c) ® Tr(e-^H^.) - ^6 Can 

imagine that left lead (L), center junction (C), and right 
lead (i?) in this model was in contact with three different 
heat baths at the inverse temperature /3l = (fcsTi) ^, 
/3c = (fcsTc)"^ and = {kBTa)'^, respectively for 
time t < to. At time t = <o, all the heat baths are re- 
moved, and coupling of the center junction with the leads 
and the interaction between the two leads are switched 
on abruptly. Now the total Hamiltonian of the lead- 
junction-lead system becomes 



-Hl + He + Hr + Hlc + 



CR 



H 



LR, 



(4) 



where Ha = ^p^Pa+^u^K°'Ua, a = L,C, R represents 
coupled harmonic oscillators, and Pa are col- 

umn vectors of transformed coordinates and correspond- 
ing conjugate momenta in region a. The superscript 
T stands for matrix transpose. Hlc = u^V^'-^uc and 



HcR = WqV'~^^ur are the usual couplings between the 
junction and the two leads, which are certainly necessary 
to establish the heat current. Now the new term repre- 
senting interaction between two leads Hlr = u^V^^ur 
will modify transmission coefficient greatly, which is our 
main interest. 

It is worth mentioning that nonlinear interaction could 
be added inside the center junction and dealt with us- 
ing self-consistent approach in the framework of NEGF, 
which has been done by many authors [3 13 • 



B. Steady state contour-ordered Green's functions 

Contour-ordered Green's functions are the central ob- 
jects in the NEGF formalism, among which the directly 
derived relation say, Dyson equation, could be readily 
transformed to all kinds of relations among the real-time 
Green's functions by Langreth theorem [201 . And many 
interesting quantities such as the current we will consider 
in the following subsection could be easily related to the 
proper real-time Green's functions. 

Steady-state contour-ordered Green's functions are de- 
fined as 



-Tr 



r{rWrc[^,"(ri)uf (r2)]}, (5) 



where p^* (s) = U (s, t^) p (to) U {Iq, s) is the steady-state 
density operator, in which time s > to introduced for 
convenience of later discussion could take any finite time 
since the switch-on time to will be let to go to — cio at 
the end in order to establish steady-state heat current. 
(ti) = U {s,Ti) u'jU (ti,s) is operator in the Heisen- 

berg picture and similarly for )■ The variables ri 

and T2 are on the contour from time s to oo and back from 
oo to time s. U{to,s) etc. are the time evolution oper- 
ators of the full Hamiltonian. Tc is the contour-ordering 
super-operator. There is a strong assumption here which 
is all we need in the whole derivation, where we assume 
steady state could be established from initial product 
state after infinite time so that all the steady-state real- 
time Green's function depend only on the difference of 
the two time arguments. This intuitively reasonable as- 
sumption is not always guaranteed and there is a specific 
example about how to establish steady-state heat current 
in Ref P|. 

After to ~oo , s t^ , and transforming to the 
interaction picture, where the total Hamiltonian Htot is 
separated into the free part Hq = Hl + He + Hr and the 
interaction part Hint = Hlc + Hcr + Hlr, wc obtain 

(n,r2) = 



-Tr|/5(-oo)rc 



(6) 



where zi" ^ (ri) = e^^o'^^u^e n.^a'^^ is operator in the in- 
teraction picture and similarly for ^ (T2) and (t'). 
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Now the variables ri and T2 are on the Keldysh contour 
[23, K from —00 to 00 and back from 00 to —00. 
The contour variables such as ri only influence the or- 
dering of the operators under T^, and ei^°'^^ has the 
same meaning as e's^"*! with real time ti. Expanding 
the exponential to perform a perturbation expansion and 
using Feynman diagrammatic technique, we can obtain 
Dyson equations for G"^ {ti,T2) , a, (3 = L,C,R such 
as Gg^ (n, rs) = Em Ik drgfi {ri,r) l^^^G^/^ (r, T2) + 

T.i,n Ik drgfi^ (^1 ' ^) Vin^'G^.I' i^, r,) , etc. All ^ these 
Dyson equations could be symbolically lumped into a 
compact matrix expression, 



G 



gVG = g + GVg, 



(7) 



where G 
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qLL qLC qLK 
qCL qCC qCR 
qRL qRC qRR 
yLC yLR 


yRC 
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9^' 
















(1/^ = V), and 



-^r4<^.(ri)<,fc(r2)] 



(8) 

a — L^G,R are equilibrium contour-ordered Green's 
functions for the free subsystems, which are easy to cal- 
culate directly. No approximation is needed here, since 
the coupling Hint is quadratic. 



C. Generalized steady-state current formula 

Certainly, heat current flowing out of the left lead in 
steady state doesn't depend on time and based on its 
definition /f = -Ty [p^" {s) dH l [t) / dt] for < > s, we 
could simply obtain 

^r--/^^^Tr[(yG<H),J 

" |^^Tr[(K,,G<,H)^J, (9) 

where {VG^ ['^Dll denotes the LL part submatrix of 
VG< [w]. Observing the structure of Tr [{VG< [w])^^], 
we note that the size of the G^ [u] making nonzero con- 
tribution to is completely determined by nonzero en- 
tries in the symmetric total coupling matrix V . So we 
don't need the full G^ [w] which is an infinite matrix 
due to the two semi-infinite leads. According to this ob- 
servation, we choose the reduced square matrix Gf^^ [uj] 
to be the corresponding submatrix of G^ [w] determined 
by the row indexes of nonzero row vectors of coupling 
matrixes V^^ , V^^, V^-^ , V^^ plus full center part row 
indexes inside the total coupling matrix V for the rows 
of G^g^ [w], the column indexes of nonzero column vec- 
tors of coupfing matrixes y'-^^ ,V^^ ^V^^, V^^ plus full 



center part column indexes inside the total coupling ma- 
trix V for the columns of G^^^ [uj]. In order to calcu- 
late the lesser Green's function Gf^^ [uj] , closed Dyson 
equation for reduced contour-ordered Green's function 
Gred (ti, T2) is ucedcd. Equation ([T]) is the starting point 

and indeed it is also true that Gred = 5red+ffredK-edGred, 

where gred is similarly defined as Gred and Vred is the sub- 
matrix of original V after crossing out all the zero col- 
umn and row vectors except for the possible zero vectors 
whose row or column indexes are the center (junction) 
ones. Actually, Gred is just the corresponding submatrix 
of G just like Vred- 

From now on, for notational simplicity, we omit the 
subscript red of all the steady-state Green's functions 
and all the coupling matrices with the understanding that 
these matrices are of finite dimensions. 

Using the Langreth theorem [20| and Fourier trans- 
forming the obtained all kinds of real-time Green's func- 
tions, we can get 



--G'' 



^-z/lFlM 

-^/i^f^[H 



G" 



where 



( sur.a \ ^ ( sur.r \ 

[9{l:r.}) - [9{LM}) 



(10) 



(11) 



gSur,a advanced surface Green's function for the left 

lead coming from the corresponding part of the advanced 
reduced Green's function g^^^ and similarly for the re- 
tarded one. This new function plays important role for 
our generalized Caroli formula and for an interface for- 
mula to be derived below. Here, fiuctuation dissipation 
theorem [uj] = fa [uj] [g^ — 5^) , a = L,G, R. has been 
used. So is (gc)^^ — (ffc)^^ = 0' which is responsible 
for the vanishing of junction temperature dependence of 
final steady-state current formula. With respect to vari- 
ous Green's functions and specific convention of Fourier 
transform, we use the same definitions as Ref. [23 |. 

Substituting the Eq. (|10p into steady current expres- 
sion we can easily obtain 

/OO 7 
^huj{fLT,[uj] + fnT2[u;]), (12) 

Where, 

Ti M = *Tr (v^^G^cl^lGI^ + F^^G^^^fiG^i) , 

(13) 

T2 H = ^Tr (f^^G^^F^GV + ^^^G^^^f ^G?^^^ 

(14) 

Again applying Langreth theorem and Fourier transform 
to the corresponding reduced one of Eq. ([7]), we could 
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get GU = Gl^V^^g'^p'^ + Gl^V^^gT''^ and GJ^ = 
Qo^^yZCga _^ Qo^^yRCgO^^ Using the relations such as 

(^cl)^ ^ ^LC (^^) ~ fi (where the superscript f 
stands for transpose conjugate) etc., we obtain 



Ti M + Tl H =Tr [G^nL^LGlj^Tn) , (15) 

In deriving it, cyclic property of the trace was used. Fol- 
lowing similar steps, we could get 

T2 M + t; H — Tr (G^^TLGl^rn] , (16) 



Due to these properties that 



Fq, , a — i, R, 



(Grl)'^ = Gip. and {G^j^^f = G^^, it is easy to show 
that Ti + = - (Ta + Tj*) . Now we define the general 
transmission coefficient 



Tg H ^Ti M + M = Tr (GJ^^FiGl^F^J . (17) 

Since current is certainly a real number, and this prop- 
erty has been kept in the whole derivation, we have 



hill 



11"*) 



According to the definitions of retarded and ad- 
vanced Green's functions in the frequency domain, 
we know G^^ [-cj] = (G'^^ [w])^and T^l^r} [-uj] = 

(y-'^{L,R} YA) ■ Together with /l [-uj] - fn [-w] = 
— [cj] + [w], steady current can be simplified fur- 
ther to the final expression 



/•OO 7 



(18) 



Thus, it is the same as expected that Landauer-like for- 
mula still apply to this general case taking lead-lead in- 
teraction into account. And this Landauer-like formula 
with the explicit general transmission coefficient expres- 
sion (fT7|) is our central result. 

Now we need to know how to calculate G£^ in order 
for specific applications. According to the corresponding 
reduced one of Eq. ([7]), we can obtain a closed equation 
foi^ GIh 



Gr 
Li?, 



(19) 

where = ((Sa"")"' " \ a ^ L,R, 

and [v'^^''-y = Y^^ -' = V^^ ^V^'^ g}.V'^^. Since 

G'rl ~ (^Li?)^' "^ow all the quantities necessary to 
obtain general transmission coefficient Tq could be ex- 
pressed in terms of retarded or advanced form of subma- 
trix of Qred and submatrix of Yredi which are both easily 
obtained. 



D. Recovering Caroli formula and deriving an 
interface formula 



First, we recover Caroli formula for transmission co- 
efficient. In this case, coupling between the two leads 



Vlr has been assumed to be 0. Thus, similar to what we 
did in subsection III C[ we could easily derived G£^ = 
5r'''^^^G-Cfl = ^r-'^^^^G^.c^^^sr'^- Together 
with G'^^ = (G£^)^, we could immediately obtain from 
formula ini) that T [w] = Tr (GJ^^rjiG^cFi). Here, we 
should remember that all the quantities inside the trace 
now arc reduced ones. However, it is still equal to ex- 
pression ([2]), in which all the quantities could be the full 
ones, taking trace operation and the reducing procedure 
for G and V into account. How to calculate T [w] and 
apply this efficient formula to specific applications has 
been stated by many authors, e.g. j25j . 

Now we try to derive an interface formula still based on 
formula (|17p . By interface we simply mean left lead and 
right lead has been connected directly and center junction 
has been removed. Mathematically, we know V'^^ = 
and V'-^^ — in this situation. Consequently, G£^ = 



gr^V^^'G^RR and G?,^ 



{GIr)' 



G^RRV^'^gr^. 



Straightforwardly, we get the transmission coefficient for- 
mula in this interface problem [26| 



(20) 



r,H =Tr GJ^^FflGJ^^Fz. 



In order to apply this formula, still we need a closed 
equation for G^^^, which could be simply obtained to be 



RR —9r + 9r ^L^RR 



(21) 



where the reduced retarded self-energy is given by E£ 

yRLgSur^ryLR^ 



III. AN ILLUSTRATIVE APPLICATION 

The illustrative example is a one-dimensional central 
ring problem, in which there is only one particle in the 
center junction connected with two semi-infinite spring 
chain leads. In this model, the interaction between the 
two nearest particles inside the two leads also exists 
taken into account as V^^. Thus, the form of the to- 
tal Hamiltonian is the same as (|4]) with A'q , a = L,R 
the semi-infinite tridiagonal spring constant matrix con- 



along the diagonal and —ujf along 



2ujt 

,2 



^0' 



LC 



sisting of 2ujf 4 
the two off-diagonals, = 

K?(f = ^red = -P^'i- where /3 is the cou- 

pling strength between two leads. The on-site potential 
term ujq is necessary in establishing the steady-state cur- 
rent dynamically [2l|. In this simple case, there is an 
analytical expression for g^Q^ [uj] , a = L, R, which is 

Cq' = -^iM' ^1 = {-n±^n^ -Aujf)/{2Lul), where 
n = {(jj + jO+)^ — 2u}f — ojq and the choice between the 
plus or minus sign depends on satisfying |Ai| < 1. And 
g^r,r _ I /ft. After all these preparations, the trans- 
mission coefficient is simply calculated by the formula 

Also, there is an alternative method to deal with this 
problem suggested by Di Ventra as we mentioned in Sec- 
tion m Essentially we repartition the total Hamiltonian 
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FIG. 1. (Color online) The transmission coefficient T [w] as a 
function of frequency for coupling between leads (a) /3 = 0.2, 
(b) /3 = 0.4, (c) 13 = 0.6, (d) /? = 0.8. The results was 
calculated directly (red solid line), and by repartitioning the 
total Hamiltonian (blue circles), ojo = 0.1 oji in all cases. 



so that interaction between leads is absent. Thus, in this 
model, the form of the total Hamiltonian is still the same 
as (m but with 



2^2 



CR 
red 



-Lol 0] 
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(22) 
(23) 
(24) 



Since now = 0, we can use either the Caroli formula 
or the general one ([T7)) to calculate the transmission 
coefficient. The results of the two methods were com- 
pared in Fig. [TJ It turns out to be that they are the 
same, which justifies the suggestion of Di Ventra from 
NEGF point of view in this example. 

Probably a much efficient way to calculate the trans- 
mission coefficient in this type of noninteracting prob- 
lem is to use the interface formula (PU)) . Frequently, the 
surface Green's functions will become complex when we 
separate the total system into two parts in order to apply 
the interface formula. However, there are some efficient 
algorithms for surface Green's functions see, for exam- 
ple, Ref. [2^. Now we will show a specific application of 
interface formula (l20l). 



IV. AN EXPLICIT INTERFACE 
TRANSMISSION FUNCTION FORMULA 

Here in this section we derive an explicit expression 
for the transmission function Tj [a;] using Eq. (j20p for the 



single interface setup, i.e. the left and right lead are 
directly connected and the center part is removed. 

Let us consider that the normalized force constant for 
left and the right leads are lJ\ and w| respectively and 
the normalized interface coupling strength is Lo\r^. Also, 
onsite potential Wg to all the atoms exists to ensure that 
the steady state could be established dynamically. This is 
a quite general scenario for a one-dimensional harmonic 
chain which is useful for the study of interface effects. So 
one of the force constant matrix say is equal to + 
Ai^ where AiiT is the semi-infinite matrix with only first 
element is nonzero Ai^Tn = w^2~'^i while K}^ is the same 
as defined in the last application. Similarly for with 
replaced by In order to obtain the explicit form, 
only inputs that are required are retarded surface Green's 
function g™''''- for both the leads. G%;^^ in Eq. ^ can 
then be easily obtained from these expressions. 

Let us calculate the surface Green's function for one of 
the leads, say the left lead. Then for the right lead it can 
be obtained just by replacing with The surface 
Green's function for a semi-infinite lead when all force 
constants are the same is given as before, i.e. g|,"g''^- Now 
for this interface case we can obtain the surface Green's 
function as follows. The retarded Green's function for 
the left lead satisfies the following equation 



[(^ + iO+)2-X^]g£=./. 



(25) 



Taking ~ + AiiT into account, and using /S.K as 
a perturbation we can write 



9l 



r,0 

9l 



(26) 



Since in this case only first atom of the left lead is con- 
nected with the first atom of the right lead, the retarded 
surface Green's fimction of the left lead is just the (1, 1)*^ 
element of g£ and we obtain 



9l 



and the self-energy for the lead is given by 



yr 



12 



^12 



(27) 



(28) 



Knowing this surface Green's function and self-energy, 
we can easily obtain Tj[lj] from Eq. (|20p which can be 
written as 



C.2c,2^t2(Ai-At)(A2-A^) 



|(a;2 - -w2/Ai)(a;f -a;^2 -wi/Az) 



^12 



(29) 

where, A2 is similarly defined as Ai with ui replaced by 
1x12 ■ It has been noted that it matches exactly with the 
result in Ref. [27[ , where this expression is obtained from 

= — "^12 then we 
= 1 for uj within the 
phonon band < < 4aj^ + uJq and outside this 
region. 



wave-scattering method. Now if uif 
have perfect transmission i.e., Tj[uj] 
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V. SUMMARY 

We examine the heat current in a lead-junction-lead 
quantum system, in which coupling between the leads 
has been taken into account. After assuming ideal steady 
state could be established from initial product state, we 
rigorously derived a general Landauer-like formula in the 
NEGF framework, from which the corresponding trans- 
mission coefficient was obtained. Based on this general 
transmission coefficient formula, Caroli formula was re- 
covered and a computationally efficient interface formula 
applicable to the case in which the total noninteracting 
Hamiltonian could be rcpartitioncd was derived. Also an 
illustrative example was given as both a verification of the 
validity of the repartitioning procedure which doesn't af- 



fect the steady current value and the clarification of the 
meaning of some quantities used in the formula such as 
V^J^ etc in different situations. Finally, wc derived an 
explicit transmission coefficient formula in a quite gen- 
eral one-dimensional interface situation based on inter- 
face formula, which turned out to be perfectly consistent 
with result obtained by wave-scattering method. 



ACKNOWLEDGMENTS 

We would like to thank Lifa Zhang and Juzar Thingna 
for insightful discussions. This work is supported in part 
by a URC grant R-144-000-257-112. 



[1] L. G. C. Rego and G. Kirczenow, Phys. Rev. Lett. 81, 
232, (1998). 

[2] M. R Blencowe, Phys. Rev. B 59, 4992 (1999). 
[3] D. Segal, A. Nitzan, and P. Hanggi, J. Chem. Phys. 119, 
6840 (2003). 

[4] A. Dhar and D. Roy, J. Stat. Phys. 125, 805 (2006). 
[5] A. Dhar, Adv. in Phys., 57, 457-537 (2008). 
[6] A.Ozpineci and S. Ciraci, Phys. Rev. B 63, 125415 
(2001). 

[7] J.-S. Wang, J. Wang, and N. Zeng, Phys. Rev. B 74, 
033408 (2006). 

[8] M. Galperin, A. Nitzan, and M. A. Ratner, Phys. Rev. 

B 75, 155312 (2007). 
[9] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 

(1992). 

[10] T. Yamamoto and K. Watanabe, Phys. Rev. Lett. 96, 
255503 (2006). 

[11] C. Caroli, R. Gombescot, P. Nozieres, D. Saint- James, J. 

Phys. G: Solid St. Phys. 4, 916 (1971). 
[12] N. Mingo, L. Yang, Phys. Rev. B 68, 245406, (2003). 
[13] W. Zhang, T. S. Fisher, and N. Mingo, Numer. Heat 

Transf. Part B, 51, 333 (2007). 
[14] S. G. Das and A. Dhar arXiv:1204.5595. 



[15] M. Di. Ventra, Electrical Transport in Nanoscale Sys- 
tems, Gambridge University Press, 2008. 

[16] R. J. Rubin and W. L. Greer, J. Math. Phys. 12, 1686 
(1971). 

[17] J.-S. Wang, B. K. Agarwalla, and H. Li, Phys. Rev. B 

84, 153412, (2011). 
[18] H. Ness, L. K. Dash, and R. W. Godby, Phys. Rev. B 

82, 085426 (2010). 
[19] Y. Xu, J.-S. Wang, W. Duan, B.-L. Gu, and B. Li, Phys. 

Rev. B 78, 224303 (2008). 
[20] H. Haug and A.-P. Jauho, Quantum Kinetics in Trans- 
port and Optics of Semiconductors, 2nd ed. (Springer, 

New York, 2008). 
[21] E. G. Guansing, H. Li, and J.-S. Wang, arXiv:1105.2233. 
[22] J. Schwinger, J. Math. Phys. 2, 407 (1961); L.V. Keldysh, 

Sov. Phys. JETP 20, 1018 (1965). 
[23] J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 

(1986). 

[24] J.-S. Wang, N. Zeng, J. Wang, and G. K. Gan, Phys. 

Rev. E 75, 061128 (2007). 
[25] J.-S. Wang, J. Wang, and J. T. Lii, Eur. Phys. J. B 62, 

381 (2008). 

[26] B. K. Agarwalla, J.-S. Wang, and B. Li, Unpublished. 
[27] L. Zhang, P. Keblinski, J.-S. Wang, and B. Li, Phys. Rev. 
B 83, 064303 (2011). 



